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^ Abstract 

This paper considers the probability density and current distributions generated by a point-like, 
isotropic source of monoenergetic charges embedded into a uniform magnetic field environment. 
Electron sources of this kind have been realized in recent photodetachment microscopy experiments. 
Unlike the total photocurrent cross section, which is largely understood, the spatial profiles of 
CN charge and current emitted by the source display an unexpected hierarchy of complex patterns, 

i 1 even though the distributions, apart from scaling, depend only on a single physical parameter. 

We examine the electron dynamics both by solving the quantum problem, i. e., finding the energy 
Green function, and from a semiclassical perspective based on the simple cyclotron orbits followed 
by the electron. Simulations suggest that the semiclassical method, which involves here interference 
between an infinite set of paths, faithfully reproduces the features observed in the quantum solution, 
even in extreme circumstances, and lends itself to an interpretation of some (though not all) of the 
rich structure exhibited in this simple problem. 
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I. INTRODUCTION 



The classical cyclotron motion of a charge q in a homogeneous magnetic field 3 is a 
simple textbook problem. Hence, it may come as a surprise that electron waves emerging 
from an isotropic point source into the magnetic field environment, which are analyzed in 
this article, should display any behavior of interest. Indeed, only a few papers have been 
devoted to the subject in the literature. 

Experimentally, such sources can be realized by the interaction of negatively charged 
ions with monochromatic laser light provided the photon energy closely matches the binding 
energy (affinity) of the excess electron (photodetachment threshold). In this regime, the De 
Broglie wavelength of the emitted electron is large compared to the size of the emiting ion, 
and the absence of a long-range interaction between the photoelectron and the remaining 
neutral atomic core means that the dynamics of the emitted electron wave can be controlled 
externally, using applied electric and magnetic fields. Experimentally, modulations of the 
overall photodetachment rate of negative ions in magnetic traps near threshold were first 
observed by Blumberg et al. and subsequently explained using perturbation theory [IH3]; 
recently, the technique has been used to observe the fine structure and Zeeman splittings 
in negative ions [I]. Similar experiments measuring the influence of an electric field on 
the photodetachment cross section [5HS] were performed following theoretical investigations 
that predicted strong modulations of the detachment rate with photon energy P-fTT]. In 
semiclassical terms, the observed behavior can be ascribed to the interference of the outgoing 
electron wave with the electron returning to the emitting ion. This intuitive trajectory- 
based picture is known as closed orbit theory [T2l [13]. The detachment cross section in 
combined electric and magnetic fields was addressed in numerous theoretical papers [T1HT9] , 
and explored in at least one experiment involving parallel fields [2D] . 

Rather than the total cross section, represented by the photocurrent, we are mainly inter- 
ested in the spatial characteristics of the electron wave emitted by the point source, which 
include both its density distribution in space, and its experimentally more accessible current 
profile. Due to the force they exert on the electron, external electric and magnetic fields 
profoundly alter the simple spherical wave pattern of a free electron source. Guided by the 
classical electron trajectories, the waves refract and fold over under the influence of the fields. 
Generally, the electron can travel along more than one classical orbit from the atomic source 
to a given destination in space; the waves associated with these multiple paths overlap and 
cause marked interference. Although approximate, semiclassical methods [2TH23] based on 
these trajectories present a valuable tool in the interpretation of features observed in the ex- 
act quantum solution of the problem. For a point source of electrons with energy E, located 
at the position r', their wave function is given by the energy Green function G(r,r'; E), a 
particular solution of the stationary Schrodinger equation in the external potentials that 
represents an outgoing wave in the vicinity of the source [2"H [25] . These Green functions can 
be found in analytical form for simple field configurations, including a homogeneous electric 
or magnetic field [2"oT - [3T)] . In an electric field environment, charges classically follow free-fall 
parabolas, and their resulting spatial distribution features a regular interference pattern that 
is conveniently explained as interference between two classical paths, akin to the double slit 
model pHl I3TH33] . This theoretical prediction was confirmed experimentally by Blondel et 
al. [3H [35] who recorded the spatial current profile of electrons stripped from a negative ion 
beam using a tunable laser in an external field, thus imaging their wave function. Because 
the interference pattern is exquisitely sensitive to the electron energy E, their technique, 
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photodetachment microscopy, has become the gold standard in the precision measurement 
of electron affinities [36H39] . More recently, a similar analysis was performed for the case of 
parallel electric and magnetic fields. The simultaneous action of both forces leads to more 
complex trajectory fields with an adjustable number of interfering paths, and correspond- 
ingly more involved interference patterns that are nevertheless in excellent agreement with 
the exact quantum solution [4U1T42] . First experimental images of the electron distribution 
observed in photodetachment in parallel fields have recently become available [43J. 



In contrast, the features of the electron wave propagating from a point source in the 
presence of a purely magnetic field have received little attention. Naively, one might expect 
a less complex pattern due to the simplicity of the applied field, but it turns out that the 
electron waves now display a richer variety of structure than in the cases mentioned above. 
To illustrate this, we display a sample plot of the electron density in Figure [T] that shows 
a variety of interference phenomena at different scales. These features have been partly 
predicted by Berry jH] who pointed out the presence of an infinite sequence of caustics, 
turning surfaces of classical motion that confine the classical trajectory fields US], in 
this problem. In the naive semiclassical model, the wave function diverges at these singular 
points [45T - T49] : this unphysical behavior can be corrected using uniform approximations 
[15| H3 |50] adapted to the specific character of the singularity. In his paper, Berry [H] 
predicted the amplification of the electron waves near the caustics, but stopped short of 
calculating the electronic wave function itself. Our aim here is to systematically examine 
the properties of the electron wave from a classical, semiclassical, and quantum mechanical 
perspective, following the leads of a former study of the parallel field configuration |42j. 
We find that much of the added complexity in the magnetic case owes to the fact that the 
electron now can travel (if at all) along an infinite number of different classical paths from 
the source to any given destination, and that the electron drifts without acceleration along 
the magnetic field direction IB. As a result, the electronic wave combines properties of open 
scattering systems and closed problems with bound motion, best exemplified by quantum 
billiards [5TH51] . although some of the features we observe, such as "backflow" of the electron 
toward the source, defy easy explanations. Another observation of our study is the extreme 
resiliency of the semiclassical method which accurately renders the quantum solution even 
under adverse conditions. 

We briefly comment on the structure of this paper. In Section [XXJ we review the classical 
cyclotron motion of a charge q in a magnetic field 3 and show that the field of trajectories 
emitted by the source, apart from simple scalings in size and time, is universal. In prepa- 
ration of the semiclassical approach, we devise a method to identify the trajectories leading 
from the source to a given destination point; as a by-product, we obtain an alternative, sim- 
pler form of the caustic structure in this system first described by Berry |Hj. Studying the 
charge density, we find that the problem under consideration is classically ill-defined: Inter- 
ference and cancellation of waves is indispensible for the continuous operation of the source 
in a magnetic field. - The quantum mechanical solution of our problem is the subject of 



Section III Apart from scaling, it depends on only a single dimensionless energy parameter 
e. We determine the energy Green function by separation of the Hamiltonian in a bound, 
"magnetic" part perpendicular to B, and free propagation along the field direction, and 
obtain a generally rapidly converging series representation of G(r,r';E) that can be inter- 
preted as a sum over discrete open and closed "channels" related to the quantized Landau 
levels in the magnetic field, but bears no formal resemblance to the classical description. 
To bridge the divide in an attempt to understand the features of the quantum solution, we 
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FIG. 1. Electron wave emitted by an isotropic point source (near top of image) in a magnetic 
field. Left panel: Trajectory field, plotted as the radial distance p (horizontal) vs. the distance in 
field direction z (vertical) (dimensionless units). Right panel: Quantum calculation of the charge 
density along a radial section for an electron energy E = 50Hujl (where ojl denotes the Larmor 
frequency of the electron), corresponding to scattering into 24 open Landau channels. Dark spots 
correspond to high density. 
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take up a semiclassical, trajectory-based approach in Section [TV] and construct approximate 
solutions for the Green function and the associated charge and current densities n(r), j(r) 
in the form of infinite series, where each term represents a possible trajectory of the elec- 
tron. In the semiclassical picture, the Landau levels arise as singularities in the spectrum 
caused by constructive interference of these paths, as pointed out already by Berry [44]. 
To correct the failure of the naive semiclassical approximation near the caustics, we also 
deploy a uniform approximation based on Airy functions [23 EE] . In Section |VJ we study 
the charge and current density distributions for different values of the energy parameter e. 
We find that the cyclotron motion confines the charge to a cylindrical tube aligned with the 
magnetic field, with the source at its center; the electron streams away from the source in 
a symmetrical pattern in both directions. As Figure [T] illustrates, the charge and current 
density profiles display remarkably rich structure in all cases except for the lowest values 
of e: An underlying finely detailed interference pattern is modulated in intensity along the 
caustics, as predicted by Berry [44J , but also contains embedded arc-like superstructure that 
superficially resembles the "quantum scars" observed in closed systems [HU E2J [57] . In the 
absence of an electric field, the density and current patterns are distinct, even though they 
often show common features. However, near the energy thresholds where a new Landau 
channel opens (odd integer values of e), we encounter strongly diverging charge and current 
density profiles, coupled with a dramatic increase in the charge density associated with a 
simple interference pattern [SB]. Furthermore, "backflow," the presence of extended regions 
in space where the electron is moving toward the source instead of flowing away, as one would 
naively expect, is prevalent in this regime. — In the concluding remarks, we summarize our 
results, and briefly explore the feasibility of experimental confirmation. 



II. CLASSICAL DYNAMICS 



We first describe the motion of a charge streaming from the source (which we will identify 
with the coordinate origin r' = 0) from a classical viewpoint, and show how to obtain the 
trajectories leading to any given destination r. For simplicity, we will align the magnetic 
field vector B = Be z with the z-axis, and keep a generic value for the particle charge q and 
mass m in our calculation. We assume that q > 0, but point out that none of our results, 
with the exception of the azimuthal angle <p(t) in Eqs. and (J8|) [59], will change for a 
negative charge if q is replaced by its absolute value \q\. (For details concerning the classical 
dynamics, see one of the authors' (A. G.) Masters thesis [60] • A formally similar analysis 
applies to the problem of electron motion in parallel fields, laid out in Ref . [4TJ [42] , where 
the purely magnetic case is obtained in the limit r\ = v Q B/£ — > oo.) 



A. Cyclotron motion 



We first find the trajectories r(i) for a particle emitted under a spherical angle (9q, </>o) 
from the source. We tackle this problem using the principle of least action, 

S(r,t;0,0) = min / dt' ( ^r 2 + qi ■ A(r)) , (1) 
r M Jo ^ ^ ' 
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using the minimal coupling Lagrangian C (r(t), r(t)) with a suitable vector potential A(r). 
We choose a gauge for A(r) that conforms to the inherent cylindrical symmetry: 



A(r) = \{Bxv) = ^(-y, x ,of , (2) 
and express C (r(i), r(i)) in cylindrical coordinates (p, z, 0): 

C(r,r) = y (p 2 + i 2 + P 2 2 ) + f P 2 • (3) 

Since the coordinates z and are cyclic, and the Hamiltonian is conservative, there are three 
conserved quantities of motion: p z = mz = const, indicates free particle motion along the 
field axis, while the momentum conjugate to is the constant canonical angular momentum 
L z of the charge: 

L z = d£/d<j) = mp 2 4>+ ^"P 2 = const. (4) 

For a trajectory traversing the origin p = 0, L z must be zero, so the trajectory undergoes 
uniform rotation in the x-y plane with an angular velocity w^: 

UL = 2^ ' (5) 

the Larmor frequency of the problem [61]. Finally, the conserved energy E, given by the 
Hamilton function of the system: 

711 ( \1 i *2 i 2 12 



H(r, p) = r • p - C = — (p 2 + z 2 + p 2 (j) 2 J = const., (6) 

is simply the kinetic energy E = \mx 2 of the charge, which therefore travels at a constant 
speed Vq. (In cartesian coordinates, the Hamilton function takes the form % = ^(p — qA) 2 , 
which is useful for the quantum description; see Appendix.) 

Solving the equations of motion for a trajectory emitted under polar angle (6*0, 4>o) yields 
the familiar helical orbit: 

pit) = — sin^ol sin(oj Lt)\ , z(t) = v tcos6 , 4>(t) = (f) — u) L t . (7) 

Note that all trajectories share the same rotational motion. It is therefore sufficient to study 
the charge dynamics in the p-z space, where the path becomes a simple sine curve whose 
frequency and amplitude depend on the polar emission angle 9q. (It is convenient to dispense 
of the absolute value in Eq. ([7]), and formally allow negative values for p.) 

We point out that the trajectory field ^ is universal in the sense that for any given set 
of parameters, the motion of the charge differs only in its range and speed. We measure 
distances p = ulp/vo, z = ulz/vq in units of the maximum cyclotron orbit diameter Vq/ul, 
and introduce a dimensionless time r = oj^t. The charge then completes a cyclotron orbit 
of radius sin^o within a time period r cyc i = n: 

p(r) = sin 9q sin r , z(r) = r cos 9 , <p[r) = 4>o — r . (8) 
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FIG. 2. Graph of the dimensionless energy function e(r) (red (gray) curve). The intersections 
with the horizontal line e = 1 (blue) correspond to the possible times of flight r. Starting from 
some initial value (here, vq = 2), we find two solutions per cyclotron interval vtt < r < 1)tt, 
implying an infinite number of trajectories. (Parameters used: p = 0.5, z = 5.) 



B. Finding trajectories 



While Eq. (|8j) details all trajectories, we are interested in finding those paths that lead 
to a given target r with dimensionless coordinates (p,z,(f>). Eliminating the emission angle 
#o in the set of equations (|8| yields implicit solutions in terms of a function e(r) of the time 
of flight r: 

= ^- + |! = I- (9) 
sin t r z 

It is easily verified that there is exactly one trajectory r(t) that connects the source with the 
destination (p, z) in a given time T if one allows the energy of the charge to vary; e(r 
then indicates the energy of this trajectory in units of the fixed source energy E = \mv l 
We therefore denote e(r) the energy function in this problem, and the condition ^ is now 
seen to pick the trajectories of the "correct" energy. 

Figure [2] shows a plot of e(r) that reveals some important features of this function: It 
is strictly positive and diverges at every integer multiple r = utc of the cyclotron period 
Tcycio, giving the graph the appearance of a succession of disjoint U-shaped curves in each 
interspersed region. The minima of these curves, located at the roots of e'(r) |62j : 

* W = - 2 fe + a, (10) 



dr \ sin 3 r r 3 

form a monotonically decreasing sequence that tends toward p 2 as the number v of the 
cyclotron interval tends toward infinity. Hence, classical trajectories do not reach beyond 
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the cylindrical region p > 1. (We note, however, that the analytic continuation of e(r) 
formally admits complex-valued solutions for the time of flight r. These lead to complex 
"ghost orbits" that lack physical reality from the classical viewpoint, but acquire meaning in 
the semiclassical description as exponentially suppressed "tunneling trajectories" that make 
the classically forbidden region accessible to waves. We will use them in our simulations in 
this paper.) 

Conversely, inside the cylinder p < 1, there will always be an infinite set of trajectories 
solving our problem. Starting from some initial cyclotron orbit uq, the graph of e(r) will 
intersect e = 1 twice in every subsequent interval v > i/q, corresponding to a pair of "fast" 
and "slow" trajectories that undergo v complete cyclotron orbits before arriving at (p,z). 
Physically, these trajectories form ever more tightly wound helices leading from source to 
destination. 

In general, the transcendental equation ^ has no analytic solutions, so we rely on nu- 
merical methods to find the times of flight r v and corresponding trajectories r u (r) of the 
charge. Our fast and reliable scheme is based on Newton's method. Since all other physical 
quantities of interest can be expressed as functions of the time of flight r, finding the roots 
in Eq. (|9| is the only instance where we resort to numerical computation. We note, however, 
that in the limit v — > oo, the times r asymptotically approach a regular sequence: 

t ~ z/7T + arcsinp and r ~ [y + 1)tt — arcsinp (v — > oo) . (11) 

This approximation is helpful when studying the convergence properties of the source prob- 
lem from the classical and semiclassical perspective. 



C. Caustics 

Although an infinitude of classical paths will connect any destination point r (with p < 1) 
to the origin, the qualitative character of the solutions generally differs among these points, 
because they do not all share the same minimum amount uq of cyclotron orbits required 
to reach the destination. For each value of uq, the set of destinations forms a separate 
manifold, and a pair of trajectories is gained or lost whenever r crosses over its boundary 
to the "neighboring" manifolds requiring z/ ± 1 cyclotron orbits. These singularities in the 
trajectory field, where the smooth mapping between trajectories r(t) and the target r breaks 
down, are known as the caustic set [HJ HH] • 

In our problem, two distinct types of caustic points are present |44j : There is an infinite 
sequence of rotationally symmetric, onion-shaped caustic surfaces that are nested inside 
each other. These surfaces are related to the cyclotron motion of the charge, and the z/th 
surface [63J represents the set of destinations where trajectories requiring v prior complete 
loops go in and out of existence. For a quantitative study, we employ the scheme for finding 
solutions (Figure [2]) introduced above. Accordingly, the z/th caustic is comprised of those 
points where the intersection of the energy functional e(r) with the line e = 1 in the z/th 
cyclotron interval is lifted, so e = 1 is tangential to the graph of e(r). Thus, the caustic is 
implicitly defined through the pair of equations: 

e(r) = 1 and |(t)=0, (12) 

where vix < t < {y + l)7r. Due to the transcendental nature of the equations ^ and (10) 
the shape of the caustic cannot be represented by a closed expression p(z). However, it is 
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FIG. 3. Caustic pattern in the p-z plane for the motion of a charge in a magnetic field. During 
their vth. cyclotron orbit, the trajectories (thin curves) trace out the z^th caustic surface (bold blue 
curves). They form an onion-like nested set. Besides these surfaces, a singularity in the trajectory 
field occurs along the symmetry axis p = (dashed red line). 

possible to reorder them, and obtain a parametrization of the caustics by the time of flight 
r: 

s / sin 3 r „, . / r 3 cos r 

P(r) = \ ~. , Kr) = ±\ — , (13) 

y sin t — t cos r V T cos r — sin r 

where r runs through the interval {y + \)k < t < {v + l)7r. Note that at the lower limit 
t = {v + ^)vr of this interval, one finds p = 1 and 5 = 0. As r increases, p shrinks and \z\ 
increases, until the caustic ends in a cusp-like structure at the symmetry axis (p = 0) at 
the locations z = ±(z/ + l)vr, the distance traveled in v + 1 complete orbits. Therefore, the 
caustics are symmetric surfaces stacked inside each other, and are joined at the circle p = 1 
in the x-y plane. (An alternative, more complicated parametrization using the emission 
angle #o has been proposed by Berry [44].) The pattern is shown in Figure [3] To illustrate 
the relation between caustics and orbits further, we include a few trajectories in the graph 
(thin line). They cross the symmetry axis v times, corresponding to v complete cyclotron 
orbits, before being "reflected" off the i4h caustic surface. Note that like the trajectory 
field, the caustic set has a universal shape, only subject to scaling. 

Besides these "fold" caustics [43], I47H49] . there is an unrelated second type of caustic 
that extends along the entire symmetry axis (p = 0) of the problem. Along the symmetry 
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axis, all trajectories emitted under a fixed polar angle 9q will converge again in a single 
point, irrespective of their azimuthal emission angle 0o- Again, this symmetry lifts the 
smooth mapping between trajectories and destination point, and leads to the formation of 
the caustic. — Finally, we emphasize that the singularities are associated with the trajectory 
field only; the propagation of a charge along a given trajectory r(t) is always smooth, even 
when it traverses a caustic point. 



D. Classical density 

Having identified the individual trajectories r v (t) that lead to the destination r, we need 
to find the relative weight carried by each of these paths in order to fully establish the 
classical and semiclassical dynamics of the charge in a uniform magnetic field. To this end, 
we examine how a bundle of "neighboring" trajectories akin to r u (t) spreads as it travels 
from the source to the target. 

Consider an infinitesimally small volume dV = pdpdzdcj) = (vo/uL) 3 pdpdz d(j> surround- 
ing the destination r. Assuming that the trajectory r„(t) arriving at r is emitted under 
angles (9q, 4>o) and has a time of flight r to the destination, the trajectories leading into the 
surrounding volume will have a spread (d6o,d<f)o) in emission angle, and dt = dr/u^ in the 
time of flight. Assuming isotropic characteristics of the source, and a total flux of J charges 
per second, the number of particles emitted within this spherical angle and time interval is: 

dN u = — - — sin 9 d9 Q d<p dr . (14) 

Since the number of particles is conserved during propagation, and the source is stationary, 
dN v is also the number of charges populating the volume dV that stem from the trajectory 
family associated with r u (t). Hence, the density n v (r) = dN u /dV due to these trajectories 
is: 

= MSU^ 0(fl o T) = M S1U0O 

Ativq p o(p,z,(p) 471^5 p 



where J = d(p,z,(j))/d(9o } (j) 0j T) is the Jacobian matrix for the transformation (J8J) relating 
the initial conditions (emission angle and time of flight) to the final position. Its determinant 
is det J = sin 2 9 t cos r + cos 2 9 sin r. Using Eq. (J8J), we substitute the emission angles via 
sin^o = pj sinr and cos# — z/t: 



n„ r 



Jul 1 



Arcv o r sin r 



p 2 cos t z 2 

■ 3 +^ 

sin r t a 



2ttvqT sin r\de/dr\ 



by comparison with the energy functional in Eq. (10). As an immediate consequence, we 
infer that the classical particle density n v {r) diverges at the location of the source (where 
r — > 0), and at both types of caustics: For sinr = 0, i. e. r = un, the charge traverses the 
symmetry axis (p = 0), while de/dr = holds along the off- axis caustic surfaces according 
to the condition ( Jl2) ). 

Outside the caustic set, the density n v {r) is finite for each individual trajectory r v (t). 
However, the total classical particle density n c i(r), obtained by summation over all paths 
leading to the destination, still diverges, as there is an infinite number of slow paths. The 
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density contribution due to such a path undergoing a large number v — > oo of cyclotron 
orbits is approximately: 



M r ) ~ T-^= " oo) . (17) 



2 1/ 



(According to Eq. (11), r ~ z/7r, and sinr w p.) Summation over all v therefore yields a 
divergent harmonic series, indicating an infinite classical density throughout the cylinder 
p < 1. This result indicates that the concept of a perfectly stationary source in a magnetic 
field is invalid from a classical perspective. Remarkably, the quantum and semiclassical 
solutions to the problem do not suffer from this defect, as we shall demonstrate now. 



III. SOLVING THE QUANTUM PROBLEM 

We now turn our attention to the quantum mechanical description of the problem, and 
study the propagation of a wave under the time evolution operator U = exp(—iHt/h), where 
"H(r, p) = 2^;[p — gA(r)] 2 is the Hamilton operator for a charge in a purely magnetic field. 
Being mainly concerned with the results, we leave a sketch of the solution to the Appendix. 



A. Quantum mechanics with sources 



One conceptual difficulty that arises in the quantum problem is that the classical tra- 
jectory field ^ describes particles spreading from the origin. Thus, "new" particles are 
continuously generated there at a constant rate. This is at odds with the conventional pic- 
ture of quantum mechanics, where the equation of continuity d t n + V • j = for probability 
density n(r,t) and current density j(r, t) implies that the number of particles is conserved. 
One workaround to the problem modifies the stationary Schrodinger equation itself, and 
introduces an inhomogeneous source term into the equation. For a point-like source of "unit 
strength," the Schrodinger equation takes the form ( |A.l ): 

[S-tt(r,p)]V(r) = S(r). (18) 



A review of the source formalism, including spatially extended sources, is contained in 
Ref.[25j 

Mathematically, Eq. ( JT8| represents the concept of a Green function, and so we call the set 
of solutions ip(r) = G(r, 0; E) the energy Green function of the charge in the magnetic field. 
(See the Appendix for a discussion.) We pick a wave function with a current field j(r) that 
radially spreads from the source, matching the isotropic emission pattern we assumed for the 
classical trajectory field in Eq. (14). (The source 8(r) has no preferred direction in space, and 
therefore invariably has s-wave characteristics. The point source formalism can be modified 
to comprise emission into p-waves and higher angular momenta [53], an extension we do 
not address in this paper.) In a more traditional interpretation, this "retarded" solution 
represents the outgoing wave part of a scattering wave function solving the conventional 
Schrodinger equation. Interestingly, the total particle current J(E) = — |9f[G(0, 0; E)] 
emitted by the point source can be directly read off the imaginary part of G(r, 0; E) at the 
location of the source, as shown in the Appendix. 
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If the motion of the charge is restricted to a two-dimensional plane, the problem of 
propagation in a magnetic field admits a closed-form quantum solution [30l 165] , while in the 
three-dimensional problem, the Green function takes the form of an infinite series [2H1 EE] , 
as outlined in the Appendix. 



B. The magnetic Green function 

In order to establish the energy Green function G(r, 0; E) for a charge in the magnetic 
field in analytic form, we note that the Hamiltonian operator %(r, p) = ^[p — <?A(r)] 2 is 
a sum of two commuting parts, a perpendicular operator %± which details the cyclotron 
motion in the magnetic field, and a longitudinal operator %\\ covering the free motion in 
field direction. (This is analogous to the classical case, where the energies in transversal 
motion E± = \mv\ sin 2 # and motion parallel to the field E\\ = \mv\ cos 2 9q are sepa- 
rately conserved.) G(r,0;E) then can be expressed as a sum over products of the various 
"transversal" eigenfunctions of a charge in the magnetic field, which are arranged in Landau 
levels Ei = (21 + 1)Hujl [67], with the simple free-particle Green function in one dimension, 
evaluated for the matching energy E\\ = E — E\. A derivation of the magnetic Green function 
using this product approach is sketched in the Appendix. (A similar procedure leads to the 
Green function for a charge in parallel electric and magnetic fields [151 113 112] , from which 
it emerges in the (non-trivial) limit 8—^-0.) For an overview of methods to calculate energy 
Green functions, we refer to Ref. [66], which contains an alternative derivation of G(r,0; E) 
in the magnetic field environment. Yet another approach is due to Gountaroulis [2"9"] . 

Unlike the classical trajectory field, the wave function G(r,0;E) is no longer universal, 
but depends on a dimensionless energy parameter e, measured in terms of the ground state 
energy hu>L of the charge in the magnetic field: 

e = E/hw L . (19) 

Alternatively, the parameter e can be interpreted as the ratio e = 2^p cyc \ /\ of the clas- 
sical cyclotron radius p cyc io = Vq/(2ujl) and the De Broglie wavelength of the charge A = 
27ih/(mv Q ); large values of e imply that the quantum length scale A is small compared to 
the classical length scale p cyc i . In terms of e and the scaled coordinates p and z, the Green 



function has the series representation (see Eq. (A. 16) in the Appendix): 




,2iy/e(e-2l-l)\z\ ~2y/e(2l+l-e)\z\ 



G(r.0:E) = \ >: L J (2e / P)_r===__- > ' L t (2ep 2 



2l+l>e 



21-1) 2 ,tt~ IV ~" V*(2Z + l-e) 



(20) 

where k = \f2mE/h is the wave number, Li(u) denotes a Laguerre polynomial, and the 
sum runs over all Landau levels (/ = 0, 1, 2, 3, . . .). As the individual terms generally drop 
exponentially once I > e/2, the series converges rapidly, and G(r,0;E) can be accurately 
evaluated numerically. The only exception is the perpendicular plane z = containing the 
source. 

Finally, we point out that the Green function G(r, 0; E) is not defined at the Landau 
levels ei = 21 + 1 itself. In the vicinity of these energies, G(r, 0; E) grows indefinitely. 
Mathematically, this relates to the observation that the Green function as the configuration 
space representation of the resolvent operator (E — H) _1 will diverge at eigenenergies E\ 
in the discrete spectrum of %. (See also the discussion in the Appendix.) Although there 
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are no bound states in the full scattering problem, the periodic motion of the charge in the 
perpendicular x-y plane caused by the magnetic field, and the ensuing quantization of the 
energy E± in transversal motion into Landau levels, becomes manifest in the singularities 
ofG(r,0;£). 



C. Quantum charge densities and currents 

From the energy Green function G(r, 0; E) we obtain both the probability density of the 
electron in the emitted wave: 

n qm (r) = |G(r,0;£)| 2 , (21) 

and the probability density current j qm ( r )' which yields the rate j- da at which charges would 
impinge on a detector area element cfa(r) (where a is the normal vector to the detector 
surface) : 

j qm (r) = ^ [G(r, 0; E)* (p - qA(r)) G(r, 0; E) + c.c] . (22) 

Unlike in the presence of an electric field S [40, 42J, where continuous acceleration tends 
to align the current profile with the density profile far from the source, the charge and 
current distributions retain their individual character in a purely magnetic environment, so 
the vector field j qm ( r ) warrants a detailed study. 

Given the symmetry of the problem, it is convenient to extract the radial component j p 
and a component j z parallel to the field: 

j z (r) = -S3 [G(r, 0; E)*d z G(r, 0; E)] , Jp (r) = -Q [G(r, 0; E)*d p G(r, 0; E)] . (23) 
m m 

Since G(r, 0; E) is a function of p and z only, the remaining azimuthal component is 
related to the vector potential term in Eq. (22), and therefore proportional to the particle 
density n qm (r): 

j (r) = -v p\G{v^-E)\ 2 , (24) 

indicating uniform rotation of the particle wave with the Larmor frequency ul, familiar from 
the classical dynamics ([7]). We will examine these quantities further in Section |VJ 

Finally, the total current J(E) can be read off the Green function. Applying relation 



(A. 3) in the Appendix, we find that only the first sum in Eq. (20), corresponding to the 



open scattering channels, contributes to the current: 



J(E) = J ircc (E) n o; ^ ' (25) 



n 2/-!- 1. . V' 6 l e 1) 



Here, Jf ree (E) = mk/(nh 3 ) is the current (A. 7) emitted by a free particle source in the 



absence of the magnetic field. Note that the total current J(E), like G(r, 0; E) itself, diverges 
at the energies e/ = 21 + 1, i. e., whenever a new scattering channel "opens" [68]. As our 



sample calculations will show, the simplicity of the result (25 ) belies the intricate and unusual 
structure present in the charge and current distributions. 
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IV. SEMICLASSICAL ANALYSIS 



In order to gain insight into the features of the quantum density and current distributions, 
we now embark on a semiclassical study of electron dynamics in the magnetic field, based on 
the helical trajectories we identified in Section [XTJ While not an exact method, we find that 
the semiclassical model faithfully reproduces the quantum results, in particular when used in 
conjunction with the uniform approximation which lifts the singularities of the semiclassical 
method that occur at the caustics. 



A. The semiclassical wave function 

To start, we assemble the semiclassical wave function ip sc (r) in an intuitive, step-by-step 
approach. (For formal reviews of the topic, see e. g. Refs. I2TH231 ) The basic idea is to 
assign a wave ip u (r) to every trajectory r u (t) that leads from the source to the destination 
r; the semiclassical wave function itself then is the sum of the individual wave amplitudes 
associated with the various paths r„(i): 

oo oo 

^(r) = j>„(r) = E hM p )l ^ ' ( 26 ) 

Note that ip sc {v) involves an infinite sum in our problem. We now proceed to define the 
modulus |?/v( r )| and the phase $^(r) of each contribution. 

The modulus ip u ( r ) is chosen so that the trajectory on its own contributes its classical 



weight n v {v) (16) to the particle density. Since |^(r)| (21) yields the probability density, 
we set |V^( r )l = \A^( r )- We note that the classical density diverges at the caustic set (the 
symmetry axis, and the nested turning surfaces displayed in Figure [3]). Therefore, ij} sc (r) 
will show unphysical behavior near the caustics. Uniform approximations (see below) are 
available to remove the divergence, and replace it with a smooth transition of the wave 
function across the caustic. 

For the phase 3> w (r), we adapt the De Broglie relation p = Kk valid for a free particle to 
the accelerated motion of the charge in the magnetic field. The classical momentum p„(r) 
then becomes position-dependent, and we postulate that the local change of the phase in 
space is again given by De Broglie's relation: V$^(r) = k„(r) = p u (r)/h. Integrating along 
the trajectory from the source to the destination yields a phase difference: 

A$dyn (r) = l _J^ p y ( r ') . dT f = * Wv ( T} o; E) , (27) 
that is proportional to the classical Hamilton- Jacobi action functional W u (r, 0; E) [69]. Using 



p = mr+gA(r), we rewrite Eq. (27) as a temporal integral, and insert the equation of motion 
([7]). We integrate, switch to dimensionless coordinates p and z, eliminate the emission angles, 
and find [70]: 

A$t yn (r) = t ( P 2 cot r + z 2 /t + r) , (28) 



where e = Efifuxii) is the quantum mechanical energy parameter (19). 

Beside this "dynamical" contribution to the phase, one must also consider a discrete 
correction that traces back to the evolution of the particle density along the trajectory. 



Eq. (15) shows that the classical density ?v(r) is a function of the determinant of the 
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Jacobian J = d(p, z, 4>)/d(9 , <p , r), where a zero of det J implies a singularity in the 
density. For the semiclassical analysis, n u (r) is an analytic function of detj', and changes 
its sign at each simple root of det J7, i. e., whenever the trajectory runs through a caustic. 
Taking the square root of n v (r) to find the modulus of the semiclassical wave function, each 
sign change in n v {r) translates into an additional factor —i, which we include into the phase 
$i, as a discrete shift of — | [7T]. The Maslov index \i v denotes the number of sign changes 
of det J along a trajectory r„(i), so the total phase accumulated from to r is: 

A$„(r) = A$^ n (r) - V . (29) 



The Maslov index \i v can be read off Eq. (16) in a straightforward manner. Whenever 



the charge has completed a cyclotron orbit and returns to the symmetry axis, sinr = 
holds, and det J changes sign. In addition, within each cyclotron interval, de/dr drops 
monotonically from +oo to — oo, implying one more simple root of det J which corresponds 
to the turning point of the path on the fold caustics. Hence, with every completed cyclotron 
orbit, the Maslov index grows by two. In the final arc of the orbit, when the trajectory 
reaches its destination, the path arrives either before touching the caustic, or afterwards, 
corresponding to the "fast" and "slow" solutions in this cyclotron interval shown in Figure |2j 
respectively. (Note that these solutions differ in the sign of de/dr.) Hence, the Maslov index 
\i v for a solution in the z/th cyclotron interval z/7r<r<(z/ + l)7ris: 

{iv ("fast" path, de/dr > 0) , 
^ ~ 1 2u + l ("slow" path, de/dr < 0) . ( ' 

Note that the semiclassical wave function ^> sc (r) can be expressed as a function of the times 
of flight r v of the trajectories. 

In order to compare the quantum result to the semiclassical approximation quantitatively, 
we finally need to fix the previously unspecified emission rate J of the classical source in 



Eq. (14). For this purpose, we identify J with the quantum mechanical current Jf ree (E) 



mk/{ixh ) (A. 7) emitted by a unit source of free particles in the absence of the magnetic 
field (see Appendix). 

To improve the semiclasssical approximation, we also include tunneling trajectories into 
our calculations. For these, conjugate complex solution pairs for the times of flight r v 
are obtained. Their semiclassical contributions follow from complex continuation of the 



expression for the dynamical phase (28); only the physically acceptable solution which leads 



to an exponentially decaying wave function is included with the sum (26) 



B. Convergence properties 



Recall that the number of classical trajectories connecting source and destination in the 
magnetic field environment is infinite, and that the concept of a stationary emitter is ill- 
defined in a purely classical description. In the semiclassical picture, the wave function 
, 0sc( r ) therefore becomes an infinite sum (26), and it is of interest to study its convergence 
as a function of the particle energy e. For this purpose, it suffices to consider the behavior 
of the individual waves ip u making up the series in the asymptotic limit v — > oo. 

In the long-time limit, the classical density n v (v) (17) drops inversely with the number of 



cyclotron orbits v. For the complex series (26), we now additionally inquire into the limiting 



15 



behavior of the phase A<& u (29) for large v. Since sinr ~ p holds in the long-time limit, we 
find cot r ~ ±a/1 — p 2 / p, so the phases for the "fast" and "slow" trajectories asymptotically 
approach the values: 

. e f + arcsin p + p\/l — p 2 ) — vtx ("fast" path) , 

A$ v ~ <! } V / v (31) 

( [y + l)ir — arcsin p — py 1 — p 2 ) — [y + ^) n ( "slow" path) , 



as v — > oo. (Note the second contribution arising from the growing Maslov index p u (30).) 
In either case, as we increment the number of cyclotron orbits v — > v + 1, the phase asymp- 
totically increases by an overall amount 7r(e — 1). Thus, the semiclassical series has the 
asymptotic form: 

^sc(r) ~ (C fast (p) + C slow (p)) ' ( 32 ) 



where the prefactor is a function of the lateral distance p only, and can be read off Eqs. (17) 



and (31) 



In the mathematical literature, the complex sum in (32) is known as a periodic zeta 
function F[^(e — 1), |] [72]. Unless the phase increase in the exponent is a multiple of 
27r, the sum is alternating, and conditionally convergent. Otherwise, the sum is real and 
divergent, growing with the square root of the summation limit. This happens whenever 
e = ei = 21 + 1 is an odd integer. Hence, the semiclasssical approximation reproduces a key 
property of the quantum solution: Unlike the classical density, ipsdr) is well-defined, unless 
the energy of the particles coincides with one of the Landau levels in the magnetic field. 



In practice, convergence of the series (32) is slow, in particular in the vicinity of the 



Landau level thresholds where e almost matches an odd integer value. In numerical simula- 
tions, we found it necessary to adjust the number of trajectories included in the summation 
to achieve good agreement with the quantum solution. We experimented with a number of 
sophisticated schemes to accelerate convergence of the semiclassical series, but found no con- 
sistent superior performance compared to a "hard cut-off" in the summation. We therefore 
adopted this simple method in our simulations below. 



C. Semiclassical density and current 



Once the semiclassical wave function \l/ sc (r) (26) is established, the analysis proceeds in 



the same vein as in the quantum case. As in Eq. (21), the approximation to the particle 
density n sc (r) = |^ sc (r)| 2 is given by the absolute square of the wave function, and now 
becomes a sum over all pairs of trajectories: 



ru 



E 

ct,/3=l 



n a [r)np[r)e 



i[* a (r)-*,)(r)] 



(33) 



It is worth noting that this sum, like the semiclassical wave function itself, is conditionally 
convergent. As a result, it cannot be reordered into a classical density n c i(r), given by the 
diagonal terms with a — f3, and "interference terms" with a ^ (3, as the classical density 
diverges (see Section In]) . 



For the current density j sc (r), we start from the quantum expression (22), and replace the 

We 



operator p — q , A(r) with its classical counterpart, the kinematic momentum mr{t) 
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obtain again a sum over pairs of trajectories, now weighted with the mean particle velocity 
v(r) at the destination point: 

_ oo 

jsc(r) = ^ E /a^[v a W+^(r)]e^«-^Wl . (34) 

a,P=l 



This is suggestive of the classical relation j = nv, but we note again that the sum (34), due to 
its conditionally convergent nature, cannot be reordered at will. In fact, the current flowing 
from a source in a magnetic field features rather counterintuitive behavior, as discussed in 
Section [V] 



D. Uniform approximation 

Because the classical density n c i(r) becomes singular there, the "primitive" semiclassical 



approximation (26) to the wave function is bound to fail near the caustics. Still, it is possible 
to find higher- level uniform approximations, based on classical trajectories, that correct the 
divergent behavior. The idea behind these approximations is that caustic points share 
generic types of divergence ("catastrophes") [45T - I49] . and therefore a wave solution valid 
in the vicinity of any such point provides a template for the solution at all related points. 
For the simplest catastrophe, the fold-type caustic, the prototype solution [45j S7J [50] is an 
Airy function Ai(u) [56], the wave solution for a quantum particle "turning around" under 
a constant force in one dimension [74] . 

Fold caustics are associated with pairs of trajectories coalescing and disappearing, 
which occurs in the magnetic field problem at the nested onion-shaped turning surfaces 



parametrized by Eq. (13), and displayed in Fig. 3 Hence, in the vicinity of these points 



the contribution of the affected pair of trajectories to the semiclassical wave function (26) 



should be replaced by an Airy function, with an argument appropriately matched to the 



change in their dynamical phase A$^ yn (27). In our simulations, we adapted a technique 



developed in Ref. 221 and assigned a combination of an Airy function and its derivative to 
each pair of "fast" and "slow" classical paths within a cyclotron period. Thus, the uniform 
approximation used is an infinite sum of Airy functions. 

The expansion in Airy function yields excellent results near the turning surfaces, but is 
not suitable for the second type of caustics encountered in our problem, the "focal line" 
p — 0. Bundles of trajectories periodically converge upon this line under all angles <p, 
and their interference yields a Bessel function of order zero Jq{u) as an amplitude profile. 
(The cylindrical cusps that form the joints between turning surfaces and focal lines need 
to be considered separately [75]. Neither focal lines nor cylindrical cusps are "generic" 
catastrophes in the mathematical sense, but are commonly encountered in systems with 
cylindrical symmetry.) Hence, the wave function at small p — > is more appropriately 
described by a superposition of an infinite number of such Bessel functions. While interesting 
in its own right, we did not attempt to model the wave function near the symmetry axis in 
this way. 

V. RESULTS 

In this section, we will present the findings from numerical studies for a representative 



set of values for the energy parameter e (19). Considering the simplicity of the setup, 
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a surprisingly rich set of features is found in these simulations, some of which defy easy 
explanation. Another objective of our studies is to assess the performance of the semiclassical 
and uniform approximations in comparison to the exact quantum results. We have pointed 
out that the semiclassical method operates here under challenging conditions, with an infinite 
number of classical paths present, wave functions represented by conditionally convergent 
series, and the classical counterpart of the problem being altogether ill-defined. However, 
we find that in all cases studied, the semiclassical technique produced reliable results, with 
an accuracy only limited by the time and depth allotted to the computations. 



A. Density profiles between Landau levels 

For our first case study, we choose e = 50 as a value for the energy parameter. This places 
the electron energy right between the 24th and the 25th Landau level, whose thresholds are 
located at 49 Hul and 51 hwi, respectively. Hence, 24 open scattering channels contribute 



to the electronic current, while the quantum wave function G(r,0;E) (20) has additional 
evanescent components that modify the density and current profile near the source. Since 
we pick a value of e removed from the thresholds 6/ = 21 + 1, the various scattering channels 
have comparable contributions, none clearly dominating the others. Hence, this example 
stands for a "regular" situation, with results "typical" for most values of e. 

We first examine the charge density n(r) generated by the source. In Figure |4j we plot 
color-coded maps of the density in the p — z plane, i. e., a cut through the three-dimensional 
distribution parallel to the magnetic field axis that contains the source. The figure shows 



the results of a semiclassical calculation (left panel) using Eq. (33), a simulation using the 



uniform approximation detailed in the previous section, based on the same trajectories (cen- 



ter panel), and the exact quantum distribution (21) (right panel). For better display, we 
plotted the "integrated" density 2irpn(p, z), summing up the contributions for various az- 
imuthal angles <ft. Dark (blue) patches in the images correspond to high density, white to low 
density. Owing to the cylidrical symmetry inherent in the problem, a section perpendicular 
to the magnetic field axis yields an interference pattern of concentric fringes that does not 
reveal any additional information. Figure 5 shows such a density distribution, where the cut 
is taken along the bottom edge of Figure 4 

The semiclassical calculations in the panels were performed with an arbitrary time-of- 
flight cutoff at r = 2507T, corresponding to at most 500 classical orbits, to keep the time of 
calculation reasonable [76] . Nevertheless, the semiclassical methods reproduce the quantum 
density distribution very well. As expected, the "primitive" semiclassical computation fails 
at the location of the caustics, where the calculated density diverges, leading to linear 
"streaks" in the image. The uniform approximation corrects this unphysical behavior, and 
provides a density map that almost perfectly matches the quantum result. Note that the 
quantum image clearly displays a strong increase in density along the location of the first 
"onion-like" caustic (see Figure |3|, and a set of interference fringes running parallel to it, 
despite the fact that the quantum calculation never invokes the concept of trajectories. This 
illustrates the power of the semiclassical method to explain features of the quantum solution. 

For a quantitative assessment of the semiclassical approximation, we also calculated a 
radial profile of the probability density, taken at z = 3.3, corresponding to the bottom edge 
of Figure |4j and representing the circular fringes in Figure [5] Here, we increased the cutoff 
time to t = 250007T, and thus added the contributions of up to 50,000 classical orbits for 
the semiclassical and uniform approximations. The resulting densities in "natural units" 
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FIG. 4. Integrated charge density 2Trpn(p, z) emitted by an isotropic point source in a magnetic 
field (vertical direction), along a radial cut in the p — z plane, for dimensionless energy e = 50. 
Left panel: Primitive semiclassical calculation (including tunneling trajectories). Center panel: 
Uniform approximation based on Airy functions. Right panel: Exact quantum result, evaluated 
using Eq. (20). Up to 500 trajectories have been summed for the approximations. The dark 



"streaks" tracing the location of the caustics in the left panel are caused by the divergence of 
the semiclassical approximation there, which otherwise agrees with the quantum simulation (right 
panel). The uniform (center) and quantum results (right) are visually indistinguishable. — The 
images cover the range < p < 1.1 in horizontal direction, and —1.1 < z < 3.3 in vertical direction. 
The source is conspicuous at the upper left edge of each image. Dark spots correspond to high 
density. All units dimensionless. 
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FIG. 5. Integrated charge density along a slice perpendicular to the magnetic field for energy 
e = 50, at a scaled dimensionless distance z = 3.3, corresponding to the bottom edge of the radial 
density maps in Figure |4j The frame extends 1.1 units from the center in x- and y-direction. Due to 
cylindrical symmetry, the quantum simulation shows circular interference fringes. For comparison, 
the concentric circles indicate the intersections with the 2nd and 3rd turning surface (inner and 
outer short dashed lines) , and the boundary of classical motion p = 1 (long dashed line) . 



[mk/(47reh 2 )} 2 derived from the free-particle Green function (see Eq. (A. 5) in the Appendix) 
are plotted in Figure [6| Both approximations quantitatively coincide with the exact result 
sufficiently far from the caustics. Whereas the primitive semiclassical method shows diver- 
gence at the intersections with the turning surfaces for p > 0, the uniform approximation is 
not affected and provides results virtually indistinguishable from the quantum calculation. 
Both methods deviate from the exact result in the vicinity of the symmetry axis p — > 0. 
This is not surprising, as p = is itself part of the caustics, and a uniform expansion into 
Airy functions is not appropriate for the focal line structure there. 

In a second step, we repeat a similar set of calculations for the current density j(r). We 
first concentrate on the current density component j 2 (r) aligned with the field direction, 
and perform again a comparison of the semiclassical and uniform approximations (34) with 
the quantum result (23). A current map, using the same parameters as the density map 
(Figure [i]), is displayed in Figure [7| We find again good agreement of the three approaches, 
except near the caustics, where the semiclassical approximation fails, as expected. While 
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FIG. 6. Radial profiles of the integrated charge density 2Trpn(p) at a vertical distance z = 3.3 
from the source, for an energy e = 50, using units of [mk/(4ireh 2 )] 2 ( |A,5 ). The plot displays the 
primitive semiclassical approximation (33) (thin red curve), uniform approximation (dashed green 



curve), and exact quantum result (21) (bold blue curve), as shown in Figure [4] (bottom edge of 
maps); the blue curve correponds directly to the density cross section displayed in Figure [5j Up 
to 5 • 10 4 trajectories have been included in the semiclassical calculation. All three methods yield 
good agreement, except that the semiclassical result diverges at the intersections with the turning 
surfaces, a failure that does not affect the uniform approximation. Both methods based on classical 
orbits are unable to reproduce the quantum density profile near the symmetry axis (inset detail 
plot). 



differing in detail, the current distribution shares the same qualitative features observed in 
the charge density distribution, including enhancement and interference along the promi- 
nent onion-type caustic. The most striking difference between the two maps is that the 
flow of particles is reversed in the upper part of the image (upward currents are encoded 
in red, downward currents in blue), simply confirming the expectation that the particles 
stream away from the source in either direction. Again, a more quantitative comparison is 
undertaken in the current profile displayed in Figure [81 corresponding to the bottom edge of 



the current maps in Figure [7| computed using the same conditions as in Figure 



(which displays the current density in "natural" units of mk 3 /(lQ7i 2 e 2 h 3 ) (A. 6), extracted 



61 The plot 



from an analysis of a free-particle source) confirms the observations about the convergence 
of the three calculations we made before. Careful scrutiny reveals an interesting detail: In 
the center of interference minima, j z {f>) drops below zero, indicating reversal of the flux 
in the upward direction. All three methods of computation agree on this counterintuitive 
"backflow" phenomenon. 
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FIG. 7. (Color online) Integrated current density in field direction 2tv pj z (p, z) in the p — z plane. 
Left panel: Primitive semiclassical calculation. Center panel: Uniform approximation based on 



Airy functions. Right panel: Exact quantum result, evaluated using Eq. (23). All parameters 



are as in Figure |4j Dark spots correspond to high current density; blue (red) indicates downward 
(upward) orientation. 
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FIG. 8. (Color online) Radial profiles of the integrated current density 2irpj z (p) for a distance 
z = 3.3 and energy e = 50, in units of mk 3 / (167T 2 e 2 h 3 ) ( |A.6 ). As in Figure [6j the primitive 



semiclassical approximation (34) (thin red curve) displays singularities at the intersection with 



caustics. The uniform approximation (dashed green curve) and exact quantum result (21) (bold 
blue curve) are in excellent agreement, except in the vicinity of the focal line p = 0. Note that all 
three methods indicate a reversal of the current at the center of some interference minima. 



Finally, we also performed quantum calculations for the radial current component j p {r) 
(23), displayed as a map in Figure [9] (left panel). Here, the current is seen to initially 
stream away from the axis p = in the vicinity of the source, but ultimately a complicated 
pattern of alternating inward flows (red) and outward flows (blue) ensues that is difficult to 
explain from a classical point of view. To illustrate the streaming pattern of the charge, we 
have combined the radial and parallel components of j(r) into a "quantum flow map" (right 
panel in Figure [9]). Here, brightness corresponds to the magnitude of the current, whereas 
the color space indicates the direction of the current probability vector in the p — z plane, 
with red indicating flux to the right, blue to the upper left, and green to the lower left. 
Transport occurs parallel to the caustic surface, but also along distinct paths of unexplained 
etiology that criss-cross the classically allowed sector from the symmetry axis (p = 0) to the 
outer limits of motion (p = 1), visible as orange and teal bands in the image. (No plot of 
the azimuthal current density j$(r) is provided, as it is simply proportional to the particle 
density n(r) (24). Classically, the whole distribution rotates uniformly with frequency w^.) 
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FIG. 9. Maps of the integrated current density vector 2irpj(p,z) in the p — z plane for energy 
e = 50. Left panel: Radial component j p (p,z) of the current density. Blue indicates outward flow 
from the symmetry axis, red inward flow; dark spots indicate high current. — Right panel: Color- 
coded map of the current flows in the p — z plane. Brightness corresponds to current intensity, hue 
to the direction of the current (red - to the right, green - to the lower left, blue - to the upper left). 
Charge transport occurs both along the caustics, and in peculiar bands (orange and teal pattern). 
(Parameters as in Figure |4j) 
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B. Density profiles near threshold 



We now slightly increase the particle energy to e = 51.01, a value just above the threshold 
of the 25th Landau level at 625 = 51. Whereas all other terms in the series solution for 
the energy Green function G(r, 0; E) (20) undergo only gradual change, the newly opened 
scattering channel adds an outsized contribution to this sum, as the relative weight of a 
given Landau level I depends inversely on the energy e — 21 — 1 available for the motion 
in field direction. The new emission mode flushes the environment with slowly drifting 
particles, and we expect that the shape of the electronic density distribution n(r) resembles 
the density profile of the dominant eigenstate ipio{f>) (A. 11) with index I = 25 of the charge 
in the magnetic field (see Appendix). The influence of the new scattering channel on the 
current distribution j(r) (22) is less obvious, as the abundance of particles in this channel 
and their slow drift velocity have opposite effects on j(r). 

From a semiclassical perspective, choosing a near-threshold energy value e implicates slow 
convergence of the series (32) for the approximate wave function ijj sc (r) (26). For such values 
of e, the phase Q v (modulo 2n) changes only slowly with each trajectory so the terms in 
the series (32) alternate over long periods, leading to large fluctuations in the partial sums. 
Hence, it is necessary to include many classical trajectories in the summation to obtain 
accurate results in the semiclassical approximation. Nevertheless, we observe that for a 
sufficiently large sample size of classical orbits, the trajectory-based methods (in particular, 
the uniform approximation) are able to reproduce the at times counterintuitive results of 
the quantum simulation. 

For comparison with the results of the previous section, we again provide maps of the 
particle density n(p, z) on a section of the p — z plane, calculated using the primitive semi- 
classical and uniform approximations (left and center panel) as well as the quantum result 
(21) (right panel in Figure 10). The caustic surfaces which structure the maps for e = 50 
(Figure [4]) are barely discernible now. Rather, the density maps now present as a sequence 
of stripes parallel to the magnetic field direction. From the quantum perspective, the uni- 
formity of the image is due to the dominance of the near-threshold scattering channel, which 
imprints its lateral profile onto the density distribution. It is remarkable that the trajectory- 
based approximations perform so well in what is essentially a rendering of a single quantum 
eigenstate in the radial direction. 

To study the density distribution quantitatively, we simulated a profile of the integrated 



radial density, taken at the bottom edge (z = 3.3) of the maps, and plotted it in Figure 11 



The peak density, near the outer classical limit, is more than thirty times larger than in the 
"regular" case e = 50 (Figure [6J , and the density profile strongly resembles the density dis- 
tribution \ipi (p)\ 2 for the Landau level I = 25 (shown in the inset) as predicted, even though 
the height of the individual maxima is modulated in the actual profile. The semiclassical 
(thin red curve) and uniform approximation (green dashed curve) again reproduce the shape 
of the quantum result (bold blue curve), save for their known failures at the caustics. 

Unlike for e = 50, the current density map j z (p, z) associated with the electron wave near 
threshold bears little resemblance to the charge density map n(p,z). We calculated j z (p,z) 
in the p — z plane using the quantum result (23) and the trajectory-based approximations 
(34) akin to Figure [7j The result, shown in Figure 12, exposes a distinctive "checkerboard" 
pattern. This pattern is partly defined by a sequence of node lines running parallel to the 
field axis which clearly trace back to the zeroes of the wave function ipio(p) for the dominant 
scattering channel (I = 25), prominent in the density plot (Figure 10). However, now the 
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FIG. 10. Integrated charge density 27rpn(p, z) in the p — z plane, for energy e = 51.01. Shown: 
Semiclassical (left panel) and uniform (center panel) approximations, and exact quantum result 
(right panel). Other parameters as in Figure g The banded appearance is due to the dominance 
of the scattering channel associated with the Landau level I = 25. 
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FIG. 11. Radial profile plots of the integrated charge density for energy e = 51.01, at a distance 



3.3 from the source (bottom edge in Figure 10), using the same units as in Figure km Thin red 



curve: semiclassical approximation; dashed green curve: uniform approximation; bold blue curve: 
quantum result. The density profile \ipio(p)\ 2 (A. 11) of the dominant scattering channel (Landau 
level I = 25) is shown in the inset for comparison. 



caustics are conspicuous in the images, and most strikingly, a rapid succession of upward 
(red) and downward (blue) current areas covers the entire plane. While the net flux of 
particles still leads away from the source, the current distribution is no longer uniform in 
direction, and almost as likely to point toward the source as away from it. While it is hard 
to pinpoint the origin of this "backflow" phenomenon, similar effects have been observed 
before in the dynamics of quantum wave packets [77], albeit not on such a massive scale as 
here. Since detectors typically are designed to absorb particles, the presence of backflow has 
profound implications for the quantum measurement process. 

For a quantitative comparison, we plot the current density profile j z (p,z) in units of 
mk 3 /(16n 2 e 2 h 3 ) at a distance z = 3.3, corresponding to the bottom edge of the maps in 
just as we did before in the regular case e = 50 (Figure pi). The result is shown 
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13 Apart from the known failures near caustics, the orbit-based approximations 



Figure 
in Figure 

again accurately trace the quantum current density. Compared with Figure |8j the magnitude 
of the current has increased about fourfold, much less than the corresponding increase in 
particle density. The plot demonstrates that backflow (negative values of the current) is a 
common occurrence, and that the reverse and forward flow are of comparable magnitude. 

To complete our comparison with the previous case e = 50, we finally also examine the 
quantum current in radial direction j p (p, z) and the current field j(p, z) associated with the 
radial and parallel components of the current for e = 51.01. The results are displayed in 
Figure [141 using the same color coding scheme as in Figure [9j In the radial current map 
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FIG. 12. Integrated current density in field direction 2npj z (p, z) in the p — z plane for energy 
e = 51.01. Semiclassical calculation (left panel), uniform approximation (center panel), quantum 



result (right panel). Parameters as in Figure 10 The conspicuous "checkerboard pattern" of 



upward (red) and downward (blue) currents indicates widespread prevalence of "backflow" toward 
the source. 
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FIG. 13. Radial profiles of the integrated current density 2irpj z {p) for energy e = 51.01; other 
parameters and units as in Figure [8j The semiclassical method (thin red curve), uniform ap- 
proximation (dashed green curve), and exact quantum result (bold blue curve) all agree in their 
prediction of sizeable ranges and magnitudes of backflow, indicated by negative current values. 



(left panel), we observe a pattern of intense, roughly horizontal bands of inward (red) and 
outward currents (blue) that almost exactly repeats in the plot of the current field as teal 
and red patterns (right panel). This means that the radial component j p (p,z) dominates 
the current distribution: The particles flow roughly back and forth from the symmetry axis, 
perpendicular to the magnetic field, with a small and variable contribution in field direction 



that gives rise to the checkerboard pattern in Figure 12 In fact, transport parallel to the 



magnetic field axis (green and purple) is almost exclusively limited to the outermost range 
of motion, p ~ 1. 



C. Exploring the classical limit 

Another avenue of inquiry is the behavior of the charge in the magnetic field if its energy 



E is large compared to the typical energy quantum hw^. Its wave function (20) then has 
contributions from many open scattering channels, and by the correspondence principle, we 
would expect that classical dynamics emerges from the quantum solution. However, since 
this problem has no meaningful classical limit, it is of interest to study the features of the 
density and current distributions obtained for large values of e. For our model simulations, 
we use e = 500, a choice for which the quantum calculation still remains easily managable. 
Figure 15 shows the integrated density distribution in this case, while Figure [T6] displays a 



map of the current density field of the charge. The color coding in the two images follows the 
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FIG. 14. Integrated current density vector 27rpj(p, z) in the p — z plane near threshold (e = 51.01). 
Left panel: The component j p (p,z) shows radial flows away (blue) or toward (red) the symmetry 
axis. Right panel: Color-coded current vector field, with brightness indicating magnitude, and hue 
representing direction (cf. Figure Note that the same underlying pattern structures both maps. 
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FIG. 15. Integrated charge density 2tt pn{p, z) in a magnetic field (oriented horizontally) along 
a radial cut in the p — z plane, for energy e = 500, evaluated using the quantum result Eq. (20). 
Canvas dimensions are < p < 1.1 and < z < 4.4, with the source located at the center of the 
left edge. Dark spots indicate high density. 



same model as in Figures [4] and |9j respectively. At this high energy, individual interference 
fringes are too densely spaced to be resolved in these images, save for the outermost range 
of classical motion [78J. Virtually all of the patterns seen are emergent "superstructure." 
Both images clearly display the first few nested caustics as curves of enhanced density and 
current, accompanied by supernumerary fringes due to the interference of the incoming and 
reflected path at the turning surface. The presence of the caustics is well understood from 
a semiclassical perspective. However, the images also show another type of superstructure, 
a multitude of thin "arcs" of various intensity that carry charge from the symmetry axis at 
the center of the distribution toward the limits of the classical range of motion and back, 
crossing through the caustics on their way. These arcs, descendants of the wide bands of 



current seen in the examples before for lower energies (Figures [9] and 14), funnel the current 
away from the source. 

The arc structure resembles in appearance the "quantum scars" frequently observed in 
billiard problems [5TJ |52j [57], which are linked to classical periodic orbits in such bound 
systems. However, the dynamics of a charge in a magnetic field is a scattering problem, 
with open trajectories, so periodic orbits are absent. Moreover, the motion is classically 
integrable, unlike the chaotic dynamics underlying quantum billiards, which renders the 
significance of the arc pattern even more mysterious. We have examined maps of varying 



energy e like those shown in Figures 15 and 16, and have found empirically that arcs are 



present at all energies, and that individual arcs persist as the energy e is increased from a 
Landau level threshold to the following threshold, heading slowly away from the source while 
undergoing fluctuations in intensity. As the energy sweeps over this interval, it appears that 
a new arc emerges closest to the source (leftmost arc in the images shown). 
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FIG. 16. Integrated current density vector field 27T/3j(/5, z) corresponding to the density map 
shown in Figure [15] Brightness corresponds to the magnitude of j, while the color hue encodes 
its direction (red - upwards, teal - downwards, yellow - to the upper left, green - to the lower 
left). Current flow occurs at the caustics, as well as along "arcs" that transport charge inward and 
outward. 

VI. CONCLUSION 

We conclude with a brief summary of our observations, and an outlook discussing their 
experimental confirmation. Notwithstanding the simplicity of the setup, the system harbors 
interesting behavior and complex features that ultimately trace back to its position at the 
junction between scattering motion common to open systems, and the periodic motion 
typical of bound systems. At the classical level, this complementarity manifests itself in 
the presence of an infinite number of orbits connecting the source to any destination inside 
the classically allowed domain p < 1, and as a consequence, an infinite classical density 
n c i(r) that renders the stationary source problem ill-defined. Thanks to interference between 
trajectories, the semiclassical model is able to lift this global singularity, and returns finite 
values for the wave function ^ sc (r) and derived quantities like the semiclassical density n sc (r) 
and current density j sc (r), except for a regular sequence of energy values E\ = {21 + l)hu>L 
(I = 0,1,2,.. .) where the divergence persists. The semiclassical results take the form of 
conditionally convergent sums over the contributions of the individual orbits. The same 
singularities occur in the energy Green function G(r, 0; E), the quantum mechanical solution 
of the problem, where they are identified as the Landau levels, the discrete eigenenergies of 
a charge confined to a plane perpendicular to the magnetic field. In general, the quantum 
solution is a rapidly converging sum over scattering channels, combinations of plane waves in 
field direction with the various radial eigenstates associated with the Landau levels. Using 
extensive simulations, we inquired into the properties of the charged particle wave emerging 
from the source, and compared the different approaches. We found that the semiclassical 
results, despite their tenuous convergence properties, were generally in excellent agreement 
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with the quantum solution, in particular when the technique was extended to incorporate 
uniform approximations based on Airy functions. The charge and current distributions 
display rich detail and a hierarchy of features, from individual interference fringes to various 
"superstructure." A nested set of shells ending in cusps on the symmetry axis is easily 
identified with the classical caustics of the system, whereas pervasive arc-like features in the 
density and current maps, reminiscent of the quantum "scars" occurring in closed chaotic 
systems, defy easy explanation. The simple energy dependence of the total current J(E) 
emitted by the source sharply contrasts with the complex spatial distribution of the current 
density j(r) which exhibits intriguing phenomena like "backflow" toward the source. 

Experimental verification of the results and simulations presented in this paper is a chal- 
lenging task, and likely not possible with state-of-the-art equipment. The established stan- 
dard for recording the current distribution of charged particle waves propagating in external 
fields is photodetachment microscopy [Ml ESI S3] • In electric field experiments, the interfer- 
ence images obtained have diameters d of several mm, and individual fringes are resolved if 
their spacing R exceeds about 100 /im [361439] . In a magnetic field, the size of the electron 
distribution is given by d = vq/ojl ^f§, and the composition of the interference pattern itself 
depends only on the dimensionless energy e = E/(fvjj£) (19). The results of Section W\ in- 
dicate that the spacing of individual fringes is approximately their ratio d/e = X/n, which 
depends solely on the De Broglie wavelength A = h/y/2mE of the charge. Hence, the energy 
should not exceed E = 2h 2 / (mR 2 ) « 1.5 • 10 -11 eV, an extraordinarily small value compared 
to the energies used in electric field experiments (10~ 5 eV). Resolution R and the radius d 
of the pattern then fix the magnetic field, B = 4h/(qRd). For d = 5 mm, the field strength 
thus obtained is minuscule, B fa 5 • 10~ 9 T. Therefore, observation of individual interference 
fringes is currently not feasible. 

Experimentally accessible values for E and B lead to fringe separations R that are much 
smaller than the resolution of the instrument. In this case, the photodetachment image 
will be governed by the superstructure imposed on the electron distribution. Because their 
location and strength shifts significantly between two Landau level thresholds, the variation 
in e = E/{huji) = d/R must be below unity in order to record the arc-like structures. 
Meeting this condition requires uncertainties in the energy distribution of the electrons, and 
drifts of the magnetic field, well below the ratio 1/e = R/d. For energies in the /xeV range, 
the uncertainties must be of order 10~ 5 , a very difficult feat to achieve. In addition, the 
unavoidable velocity spread of the negative ions in the beam [79] causes varying motional 
electric fields £ mo t that must be kept small in comparison to the Lorentz force on the electron, 
£ mo t <C VqB. Thus, unless conditions are very precisely controlled, the spread in e will wash 
out the modulations in the image due to the arc-like features. In comparison, the location 
of the caustic surfaces is rather insensitive to small changes in the parameters, and they 
should therefore be most easily recognized in photodetachment images as bands of increased 
intensity, as discussed earlier by Berry [44] . 
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Appendix: Energy Green Function of a Charge in a Magnetic Field 

In this appendix, we briefly discuss the relationship between the energy Green function 
G(r, 0; E) and the current emitted by a monochromatic point source located at the origin, 
establish the particle density and currents associated with a free-particle source as useful 
quantities for scaling the results in a magnetic field, and outline a method to find the Green 
function in a homogeneous magnetic field itself. 



1. Green Function and Current 

The energy Green function G(r, 0; E) associated with a stationary Hamiltonian "H(r, p) 
is a solution of the inhomogeneous Schrodinger equation with energy E: 

[E-n(r,p)]G(r,0;E) = 5(r) . (A.l) 

It can be shown [SO] that the Green function is uniquely defined only if E is not part of the 
energy spectrum of H. If E is one of the discrete eigenenergies of the Hamiltonian, G(r, 0; E) 
does not exist at all, whereas an entire space of solutions is available if E is a member of the 
continuous spectrum of H. This situation occurs in scattering problems without a confining 
potential, where the particles are free to leave the system. For our magnetic field problem, 
an interesting "mixed" case arises where the motion is bound in the x-y plane, leading to 
discrete Landau energy levels associated with the magnetic field, but particles are traveling 
freely along the z-axis, with a continuous range of energies. 



The degeneracy in the solutions for Eq. (A.l) in the continuous spectrum corresponds 



physically to the ability to impose boundary conditions on scattering wave functions. Since 
we study a source of particles, we are interested in the retarded Green function, the particular 
outgoing wave solution that carries particles away from the source. In the vicinity of the 
origin, this wave invariably gains the simple isotropic characteristics of a spherical s-wave 
[80] . and we therefore identify it with the radially spreading trajectory field employed in the 
semiclassical study. 



The usual expressions for the particle density n qm (r) (21) and the current density j qm (r 



(22) hold in the outgoing electron wave. To find the total current J(E) emitted by the 



source, we note that the solutions to Eq. (A.l) obey a modified equation of continuity that 
includes a source term at the origin [21)] : 



divj qm (r) = -^S [G(v, 0; E)} 5(v) . (A.2) 

While the Green function G(r, 0; E) is itself divergent near the source, its imaginary part has 
a well-defined limit which directly yields the total current J(E). We apply Gauss' theorem 



to Eq. (A.2), and integrate over a surface enclosing the origin to find: 

J(E) = -|(limS[G(r,0;£)]) . (A3) 



2. Free-Particle Green Function 



In the absence of electric and magnetic fields, i. e., for the free-particle Hamiltonian 
%free = p 2 /(2m), the energy Green function becomes an outgoing spherical s-wave. A 
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detailed analysis yields [80J: 



G free (r,0;E) = - — , (A.4) 

2nn z r 



where k = mvo/h is the wave number of the electron. For comparison, we state the particle 
density nf ree (r) and current density jf ree (^) for the free-particle source, using the cyclotron 
length units introduced in Section [II] We set f = ult/vq and note that the dimensionless 
quantity kr = 2er becomes a function of the dimensionless energy e = E/(?vjJl). Then, 

m 2 k 2 1 

nfree r = 2 2 4 , A.5 
lbir 2 e z ri i r 2 

while the magnitude of the current density is, 

J fr ee(r) = n bee (r)v = ^^3 ^ ■ (A.6) 



Finally, applying Eq. (A. 3) to the Green function yields the total current Jf me (E): 

_ , . TTik , . , . 

Jfree(E) = — . A. 7 

The characteristic growth of the source efficiency with the square root of the energy is known 
as Wigner's law 



3. Finding the Magnetic Green Function 

Finally, we briefly outline a derivation of the energy Green function in a uniform magnetic 
field. A detailed discussion, including an alternative approach to obtain G(r,0; E), is found 
in Ref . M 

We first note that the Hamiltonian H is separable into two commuting parts, a two- 
dimensional operator %± that describes the dynamics of the electron in the plane perpen- 
dicular to B, and "Hy, a one-dimensional free-particle Hamiltonian that describes the drift 
along the magnetic field axis: 

n = ^-(p-gA) 2 = n±+n\\ , (a.8) 

where H\\ = p 2 /(2m). H± is closely related to the two-dimensional quantum harmonic 
oscillator: 

n ± = ± (pi + pi) - u L L z + ^ (x 2 + y 2 ) . (A.9) 

(L z denotes the angular momentum component in ^-direction.) Since the electron motion is 
bound in lateral direction, the spectrum of 7i± is discrete, and given by the Landau levels: 

E ¥ = (21 + l)hcu L , (/ = 0,1,2,...) , (A.10) 

i. e., ei^ = 21 + 1. Each Landau level has infinite degeneracy, as the magnetic quantum 
number of the electron \i = —I, —1 + 1,-1 + 2,... can take any integer value greater or 
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equal to —I. However, only the eigenstates with p = are of interest here, since all 
others vanish at the source location. In polar coordinates, the corresponding normalized 
eigenfunctions ^o(p) read, using again the dimensionless length p = ulp/v : 



MP) = J^L^ep^e-^ , (A.ll) 

where Li(u) denotes a Laguerre polynomial [67]. Since Li(0) = 1, = l muj L/ '(Tih)} 1 / 2 S^o 

holds at the origin p = 0. Therefore, the completeness relation for the eigenstates of / H± 
can be stated as: 



OG 



£ M*> «)* = ^ E L ^ 2e ^ 2 ) = WW • (A.i2) 

Physically, the separability of the problem implies that the wave function of the electron 
is a superposition of scattering waves that occupy different Landau levels ipio, with the 
"remaining" energy En = E — Eiq associated with the drift motion in ^-direction. Indeed, 
if we introduce the one- dimensional free-particle Green function Giv(z, 0; E) as the solution 
of [E - G m (z, 0; E) = S(z), we find, 

[E-U}Mp)Gw(z,0;E-E io ) = ^o(/0 [E-E m -Un] G w (z, 0; E - E l0 ) = Mp)5(z) ■ 

(A.13) 



In view of Eq. (A. 12), the full Green function G(r, 0; E) has the series representation: 

oo 

G(r,0;E) = J2Mp)M0)*Gm(z,0;E - E l0 ) . (A.14) 



1=0 

The remaining task is to find an expression for Gm(z, 0; En). As a free-particle Green 
function, it must have the form of an outgoing wave in either direction of the z-axis (or in 
the case En < 0, an evanescent wave). A simple analysis |66j shows that: 

m e »fe|||z| > q\ 

G w (z 7 0;En) = { . . . " ' (A.15) 

-wk e l|N (^ii <o) , 



where k\\ = y/2m\E^\/h. We introduce again dimensionless coordinates, combine Eqs. (A. 12) 



(A.15), and finally obtain the full Green function: 



k ( 2iy/e(e-2l-l)\z\ -2 x /e(2i+l-e)|i| 

2 ^ 2 W<e P V^-2/-l) ^ ' V(2Z + l-c) 



(A.16) 

The summation runs over all / = 0, 1, 2, . . .. If the energy coincides with one of the Landau 
levels ei = 21 + 1, the Green function diverges. 
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